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^ ' Binary systems can be powerful sources of non-thermal emission from radio 

to gamma rays. When the latter are detected, then these objects arc known 
as gamma-ray binaries. In this work, we explore, in the context of gamma-ray 
binaries, different acceleration processes to estimate their efficiency: Fermi I, 
^ ' Fermi II, shear acceleration, the converter mechanism, and magnetic reconnec- 

' tion. We find that Fermi I acceleration in a mildly relativistic shock can provide, 

^ although marginally, the multi-10 TeV particles required to explain observa- 

■ tions. Shear acceleration may be a complementary mechanism, giving particles 

the final boost to reach such a high energies. Fermi II acceleration may be too 
slow to account for the observed very high energy photons, but may be suitable 
to explain extended low-energy emission. The converter mechanism seems to 
require rather high Lorentz factors but cannot be discarded a priori. Standard 
relativistic shock acceleration requires a highly turbulent, weakly magnetized 
downstream medium; magnetic reconnection, by itself possibly insufficient to 
reach very high energies, could perhaps facilitate such a conditions. Further 
theoretical developments, and a better source characterization, are needed to 
' pinpoint the dominant acceleration mechanism, which need not be one and the 

' same in all sources. 

_ 
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1. Introduction 

During the last decade, binary systems have turned out to be a new class of 
gamma-ray sources whose numbers are growing with the increasing sensi- 
tivity of the new instrumentation (see Refs. 1-13). Different types of objects 
pertain to this class, like microquasars, binaries hosting a non-accreting pul- 
sar, massive star binaries, and even symbiotic stars. All these sources share 
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the characteristic of hosting powerful outflows that interact with them- 
selves and their environment, and dissipate their energy partially in the 
form of non-thermal particles. Given the typical compactness of the emitter, 
the dynamical and radiation timescales are short, yielding rapidly variable 
emission that can reach high luminosities, and also very high energies. Par- 
ticularly interesting in this regard are those gamma-ray binaries that reach 
energies TeV. In some cases, like LS 5039,^*'^^ the emitting particles 
could be as energetic as ^ 100 TcV, which poses serious constraints on any 
particle acceleration model that aims at explaining the observed radiation, 
as noted by Khangulyan et al. (15). Rieger et al. (16) discussed different dif- 
fusive acceleration processes for microquasars, and also concluded that any 
mechanism responsible of the very high-energy emission should run at its 
limits. In this work, we carry out a semi-quantitative analysis of the require- 
ments and efficiency of diffusive acceleration processes (Fermi I, Fermi II 
and shear accclcration^^"^*^), the converter mechanism and magnetic recon- 
nection (e.g. Refs. 20,21), for gamma-ray binaries in general. Our goal is to 
take a first look at the problem of extreme acceleration in compact galac- 
tic sources, in which dense radiation fields, together with highly supersonic 
sometimes relativistic bulk motion, shear layers, turbulence, and magnetic 
fields, are expected. 

In Figure 1, a sketch of the general binary scenario discussed here is 
presented, showing the relevant elements that could play a role in the pro- 
duction of very energetic particles. An interaction structure is formed due 
to the presence of, for instance, a powerful relativistic outflow from a com- 
pact object (e.g. a jet or a pulsar wind) and a strong stellar wind. Two 
stellar winds could also form a similar though non-rclativistic structure. 
For the case of a jet, it will be more collimated, but jet disruption may 
also lead to a broadening of the interaction region. Powerful shocks are ex- 
pected to form at the colliding region: a termination shock in the pulsar, 
and an asymmetric re-coUimation shock when a jet is present. In the jet 
scenario, internal shocks can also occur. The contact discontinuity between 
the different flows involved is subject to Rayleigh- Taylor and/or Kelvin- 
Helmholtz instabilities (neglecting the role of the magnetic field), which 
will trigger turbulence downstream the flow, as well as mixing from the two 
different media. Despite of its complexity, the picture can be approximately 
analyzed (see, e.g., Refs. 22-24), and different acceleration site and mecha- 
nism candidates can be proposed. In the presence of strong shocks, diffusive 
acceleration could be the dominant mechanism. For magnetized, highly tur- 
bulent and diluted media Fermi II could be at work, and if strong velocity 
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gradients are present, shear may occur as well. For highly relativistic and 
radiation or matter dense environments the converter mechanism could play 
a role, whereas magnetic reconnection could take place under the presence 
of strong irregular magnetic fields. The modeling of observations tends to 
favor leptonic models (e.g. Refs. 15,25-29), although hadronic models can- 
not be discarded (e.g. Ref. 30; see also Ref. 31 and references therein). 
We will focus here mainly on electron acceleration, strongly affected at the 
highest energies by synchrotron losses, but some of our conclusions apply 
to protons as well. 




Accreting system 



Fig. 1. Sketch of the generic binary scenario. 



2. Diffusive acceleration mechanisms 

Diffusive shock or Fermi I acceleration takes place through repeated particle 
bouncing upstream-downstream of a shock front. The particle deflection is 
mediated by magnetic field irregularities of the background plasma. In each 
cycle, the particle energy gain is /S.E/E oc (us/c), where is the shock ve- 
locity. For a mildly relativistic shock and in the Bohm diffusion limit, the ac- 
celeration timescale is iacc ^ 27r {c/vsY E/qBc k. 3 (0.5 c/wg)^ -BtcV Bq^ s, 
where -EtcV is the particle energy in TeV and Bq the magnetic field 
in Gauss; Vg ^ 0.5 c would be the validity limit of the non-relativistic 
assumption. For reasonable radiation and matter densities, in a finite 
size homogeneous accelerator, dominant particle/energy losses are diffu- 
sion escape, and adiabatic and synchrotron cooling, with typical timescales 
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idiff = R^/2D = 15000 i??2 5G£^^cV s, tad - R/v - 100i?i2/WlO s, 

and = l/a^B'^E « 2,%QBq^E-1^ s, where R = Rn cm and 
V = lO^^uio cm are typical lengths and flow velocities in the re- 
gion, and D the Bohm diffusion coefficient. This yields a maximum en- 
ergy for Fermi I acceleration of i?max = i^in[i?m£{i ^mlxi -^maxli where 
Etl, - 74i?i2BG(0.5c/z;,)-i TcV, S^i, ~ mB^R^^v^^ (0.5 c/i;^)-' 
TeV, and E""^^^ ~ 12 B^^/^ (0.5 c/l's)"^ TeV. These simple estimates al- 
ready show that for slow shocks (e.g. between two massive star winds), or 
outside the range B ^ 0.1 — 1 G, it is very hard to accelerate particles 
up to > 10 TeV. It is worth noting that particle acceleration in highly 
relativistic shocks could work,^^ but requires rather speciflc magnetic fleld 
conditions downstream, such as high turbulence and relatively small mag- 
netization. This might be hard to realize in pulsar winds, in which a 
significant toroidal B-fleld is usually expected to remain downstream the 
termination shock (see however below). 

In stochastic particle acceleration or Fermi II, the situation is worse 
than in Fermi I. This process is driven only by stochastic collisions with ran- 
domly moving, magnetic field irregularities. It is a second order process, i.e. 
A£J oc (wa/c)^, where wa is the Alfven speed. Since tacc ~ (c/i^a)^ E/qBc, it 
is required that turbulent energy will be a significant fraction of the plasma 
energy and that ua c. It means that B should be around equipartition 
with the turbulence rest mass energy density. Under such a condition, if 
B ^ 1 G (i.e. optimal for vj^ ^ 0.5 c), then n ~ 100 cm~^, hardly feasi- 
ble for a compact binary. Downstream of a pulsar wind termination shock, 
relativistic Alfvenic speed may be achieved, although size and turbulence- 
energy requirements favor the region behind the pulsar with respect to the 
star. This however requires negligible stellar wind contamination and adia- 
batic losses and may be unrealistic (see Ref. 23), so a proper assessment of 
the Fermi II efliciency here requires a detailed study. In general, stochastic 
acceleration seems to be more suitable to explain extended and low-energy 
emission, in regions downstream shocks or rich in instabilities. 

Shear acceleration is, like Fermi II, a stochastic process, but relies 
on an additional global velocity gradient in the flow ~ Au/AR. In the 
mildly relativistic case, tacc ^ 3(Ai?)2/rgC = 300 Bq ARf^ E:^^y (AR = 
10^^ Ai?ii cm; rg is the particle gyroradius), and therefore the acceleration 
timescale has the same dependence on E as synchrotron. In order to operate 
then, shear has to overcome synchrotron cooling, implying B < AR^^^"^ G. 
Shear also requires of an injection process, since otherwise adiabatic or ad- 
vection (escape) losses will block it, i.e. E > 3 ARu vio Bc{AR/R)o,i TcV, 
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where (Ai?/i?)o.i means AR = R/10. In principle, Fermi I could be a good 
injection candidate, since it does no require high B to operate. Fermi II oth- 
erwise needs higher B that may render synchrotron cooling dominant over 
shear acceleration. The shear maximum energy (in the mildly relativistic 
case) is limited by = AR, i.e., r^qB ARf^ 30 Bq ARn TeV. 

3. The converter mechanism and magnetic reconnection 

In conventional relativistic (F 1) shock scenarios, charged particles are 
quickly overtaken by the shock so that they do not have enough time to 
isotropise in the upstream region. Thus, when caught up by the shock, the 
shock normal/particle motion angle will be ~ 2/r, and the energy gain, 
AE/E F^ ^^/2, will be reduced down to ~ 2. However, particles may 
have time to isotropise upstream if they managed to switch to a neutral 
state and propagate far from the shock without deflections in the S-field. 
For instance, an electron of energy E cooling via Klein-Nishina (KN) inverse 
Compton (IC) can transfer most of its energy to a scattered photon. For pair 
creation mean free paths Z ^ F^ rg, the gamma ray will pair create with an 
ambient photon far enough from the shock to allow the subsequent electron 
(or positron) to get deflected by ^ ~ 1, i.e., AE/E ~ F^/2. This effectively 
implies tacc ^ fg/c, the electrodynamical limit (e.g. Rcf. 34). If otherwise I 
is too small, then — > 2/r, i.e., the standard case. Once started, the mech- 
anism proceeds efBciently and yields a rather hard electron spectrum until 
I, roughly oc E, becomes > Fi?. After that, the electron spectrum becomes 
softer. In binary systems with bright UV stars, the converter mechanism 
could yield hard electron spectra up to ~ TeV for F ~ 100 (see also Ref. 
35). This might be the case in microquasar e^-jets before suffering mass- 
loading (caveat: gamma-ray absorption), or at the reaccelerated shocked 
pulsar wind.^^ The energy is limited by Tg ~ Fi?, although synchrotron 
losses can reduce this limit. For a pulsar termination shock in a UV stellar 
field the process cannot work efficiently, since for ^ ~ 1, a distance ahead 
the shock of F^ rg w 3 x lO^i (F/IO^)^ EtcV B^^ cm would be required. As 
shear acceleration, the converter mechanism requires an injection mecha- 
nism. In the leptonic case, E should be enough to pair-create in the ambient 
photon field (~ 30 GeV for stellar photons). 

Magnetic reconnection is perhaps the less well characterized process 
among those discussed here. Numerical calculations show that, beside bulk 
acceleration in the current sheet up to wa, non-thermal particles can be also 
produced. ^"'^ Potentially, the mechanism is fast with tacc ^ fg/c (assuming 
e ~ U, where e is the current sheet electric field), and particles may reach 
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energies limited by rg ^ AR, where AR is the current sheet size. However, 
unless the current sheet occupies a significant fraction of the source, the 
process will not yield the required multi-10 TeV particle energies. A possi- 
bility could be many reconnection sites, which could be equivalent to the 
Fermi II acceleration process. An interesting role of magnetic reconnection 
may be the dissipation of an alternating polarity S-field in the jet base (e.g., 
Ref. 37), or at the pulsar wind termination (e.g., Ref. 38), thereby possibly 
allowing further acceleration via a Fermi I-type mechanism in relativistic 
shocks (e.g., Ref. 39). 

4. Conclusions 

The present work indicates that for typical conditions expected in gamma- 
ray binaries, the production of photons with energies > 10 TeV indeed 
requires very efficient acceleration with tacc < 10 — lOOrg/c. Although less 
strictly, this conclusion also applies to protons. All this implies strong turbu- 
lence, and relatively weak magnetic fields (for electrons), as well as at least 
mildly relativistic speeds. Fermi I acceleration, although marginal, seems to 
be the best candidate in mildly relativistic outflows, but for highly relativis- 
tic flows the situation is less clear. Shear acceleration and Fermi II could 
also operate, but the latter is unlikely to help at TeV energies. The converter 
mechanism and magnetic reconnection cannot be discarded, although they 
require quite specific conditions. A better source characterization is needed 
for a proper assessment of the feasibility of all these mechanisms. 
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